THE CRITICAL TEMPERATURE FOR THE ISING MODEL ON 
PLANAR DOUBLY PERIODIC GRAPHS 



c/3 



DAVID CIMASONI AND HUGO DUMINIL-COPIN 



Abstract. We provide a simple characterization of the critical temperature for the Ising 
model on an arbitrary planar doubly periodic weighted graph. More precisely, the critical 
inverse temperature /3 for a graph G with coupling constants {Je) e eE(G) is obtained as the 
unique solution of a linear equation in the variables (tanh(/3J e )) e gB(G)- This is achieved by 
studying the high-temperature expansion of the model using Kac-Ward matrices. 



1. Introduction 



1.1. Motivation. The Ising model is probably one of the most famous models in statistical 
physics. It was introduced by Lenz in [26] as an attempt to understand Curie's temperature 
for ferromagnets. It can be defined as follows. Let G be a finite graph with vertex set V(G) 
■ and edge set E(G). A spin configuration on G is an element a of { — 1, +1} V ( G \ Given a 

positive edge weight system J = {J e ) e &E(G) on G, the energy of such a spin configuration a 
is defined by 

J4?(a) = — J e a u a v . 

e={u,v}eE(G) 

Fixing an inverse temperature (3 > determines a probability measure on the set 0(G) of 
spin configurations by 

k>A°) = "zJJgY 1 

where the normalization constant 

o-en(G) 

is called the partition function of the Ising model on G with coupling constants J. 

In this article, we focus on planar locally-finite doubly periodic weighted graphs (^,J), 
i.e. weighted graphs which are invariant under the action of some lattice A ~ Z © Z. In such 
case, 5^/A =: G is a finite graph embedded in the torus T 2 = R 2 /A. By convention & will 
always denote a doubly periodic graph embedded in the plane, while G will denote a graph 
embedded in the torus. Ising probability measures can be constructed on as limits of finite 
volume probability measures [29]. In particular, the Ising measure at inverse temperature /3 
on 5f with + boundary conditions will be denoted by ^ j3 - 

We further assume that 5f (or equivalently G) is non-degenerate, i.e. that the complement of 
the edges is the union of topological discs. A Peierls argument [32] and the GKS inequality [15| 
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20j classically imply that the Ising model on 5f exhibits a phase transition at some critical 
inverse temperature f3 c £ (0,oo): 

• for < (3 C , fj£j(a v ) = for any v £ V(Sf), 

• for /3 > /3 C , $ p{v v ) > for any v G F(^). 

This article provides a computation of the critical inverse temperature for arbitrary non- 
degenerate doubly periodic weighted graphs (Sf , J) as a solution of an algebraic equation in 
the variables x e = tanh(/3J e ). 

1.2. High-temperature expansion of the Ising model. The result will be best stated 
in terms of the high-temperature expansion of the Ising model, which we present briefly now. 
As observed by van der Waerden [35] , the identity 

exp({3J e a u a v ) = cosh(/3J e )(l + tanh((3J e )a u a v ) 

allows to express the partition function as 

Z l( G ) = ( II cosh (M)) II (1 + tanh(/3 J e )a u a„) 

eeE(G) cren(G) e=[uv]eE(G) 

(1) =( n cosh(/3j e ))2i y ( G )i n tanh (W' 

where ${G) denotes the set of even subgraphs of G, that is, the set of subgraphs 7 of G such 
that every vertex of G is adjacent to an even number of edges of 7. As a consequence, the 
Ising partition function is proportional to Z(G,x) := J2ye<ff (G) x (t)> where ^(7) = n e e7 :Ee 
and x e = tanh(/3J e ). This is called the high-temperature expansion of the partition function. 

1.3. Statement of the result. Let (&,J) be a planar non-degenerate locally-finite doubly 
periodic weighted graph. Recall that G = &/A is naturally embedded in the torus. Let Sq(G) 
denote the set of even subgraphs of G that wind around each of the two directions of the 
torus an even number of times. Set $\(G) = ${G) \ Sq{G). 

Theorem 1.1. The critical inverse temperature f3 c for the Ising model on the weighted graph 
(^, J) is the unique solution < /3 < 00 to the equation 

(2) £ x( 7 )= £ x( 7 ), 

where ^(7) = nee7 Xe anc ^ Xe = tanh(/3J e ). 

Note that the sign of X^ye<?o(G) x (^) — S 7 e,r 1 (G) indicates if the system is in the dis- 
ordered phase (when the quantity is positive) or ordered phase (when it is negative). 

Let us briefly summarize the strategy of the proof. Our main tool is the so-called Kac- 
Ward matrix associated to the weighted graph (G, J) and to a pair of non-vanishing complex 
numbers (z, w), see Definition Q] below. We show that the free energy per fundamental domain 
can be expressed in terms of the Kac-Ward determinants (Lemma I4.3j) . Using exponential 
decay of the spin correlations in the disordered phase [TJ together with a duality argument 
(Theorem 14. 4p . the free energy is then shown to be twice differentiable in each variable J e 
at any f3 ^ f3 c . The proof is completed by observing that Equation ([2]) is equivalent to 
the vanishing of the Kac-Ward determinant at (z,w) = (1, 1), which translates into the free 
energy not being twice differentiable in some variable J e . 
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Figure 1 . The graphs of Examples [HOE] and 

The main technical difficulty is to make sure that these Kac-Ward determinants never 
vanish for z and w of modulus f, except at (z,w) = (1, 1). This is achieved by showing that 
these determinants are proportional to the Kasteleyn determinants of an associated bipartite 
graph (Theorem I3.ip . a graph first considered by Fan-Wu [H] (see also [UJ). One can then 
harness [21] to show that the associated spectral curve is a special Harnack curve, and the 
desired statement follows (Lemma 14. 2p . 

This connection between the Ising model and dimers on a bipartite graph is interesting on 
its own. It was used recently by Dubedat to show the exact bozonisation of the Ising spin- 
correlations We give some other applications in the body of the text, and refer to [8] for 
further ones. Let us also mention that Kac-Ward matrices are related to s-holomorphicity and 
fermionic observables introduced in [7J (see [12] for an overview on the subject). Therefore, 
this identification of the critical inverse temperature using the Kac-Ward matrices opens new 
grounds in the understanding of universality and conformal invariance for the Ising model on 
arbitrary doubly periodic graphs. 

We would like to highlight the fact that in the case of doubly periodic coupling constants 
on @ = Z 2 , Theorem 11.11 was proved independently by Li [27] (although the statement is 
formulated quite differently). The proof in [27J relies on the paper [28] of the same author. 
The main ingredient is a sequence of transformations mapping the Ising model on Z 2 to the 
dimer model on the associated Fisher graph, and ultimately to a dimer model on a bipartite 
graph using [11] . In our strategy, the Kac-Ward matrices allow us to work on the original 
graph almost throughout, avoiding in particular the use of the Fisher correspondence. 

1.4. A corollary. A consequence of the proof is the following corollary, which identifies the 
critical point as being the unique singularity of the free energy per fundamental domain, see 
([6]) for the definition. 

Corollary 1.2. The free energy per fundamental domain is analytic in (3 for any (3 ^ j3 c . 

1.5. Some examples. We now illustrate Theorem 11.11 with several examples. 

Example 1. Consider the square lattice with coupling constants J e = 1. It can be represented 
by the toric graph illustrated to the left of Figure [U The equation of Theorem 11.11 then reads 
1 = 2x + x 2 , where x = tanh(/3). This leads to the well-known value 

P c = tanh- 1 (v / 2- 1) = -log(V2 + l). 

This value was first predicted in [23] and proved in [31] . Several alternative derivations have 
been presented, see e.g. [HEH1]. 

Example 2. Consider now the hexagonal lattice with coupling constants J e = 1 (see Figure [T]). 
This time, the equation reads 1 = 3x 2 , leading to 

p c = tanh-^v 7 ^) = ~ log(2 + ^3). 
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Figure 2. An edge e of an isoradial graph, the associated rhombus, and the 
half-rhombus angle 9 e . 

Example 3. Similarly, the equation associated to the homogeneous triangular lattice as rep- 
resented in Figure Q] is 1 + x 3 = 3x + 3x 2 . This leads to 

/3 C = tanh~ 1 (2 - VS) = - log(V3). 

Example 4. More generally, consider a doubly periodic graph ^ isoradially embedded in the 
plane, i.e. embedded in such a way that each face is inscribed in a circle of radius one, with 
the circumcenter in the closure of the face. To each edge e, associate the coupling constant 

T 1 /l + sin# e \ 

Je = ~ log ^ ' 

2 \ COs6 e J 

where 9 e 6 (0,7r/2) is the half-rhombus angle associated to the edge e, as illustrated in 
Figure [2J 

It follows from [TU1 Theorem 4.7] (see also [3]) that Yl-ye^b(G) ~ S 7 g^i(G) x (^) van ishes 
for /3 = 1. By Theorem 11.14 this is the critical temperature for the isoradial graph . The 
cases with constant angles 9 e = vr/4, ir/3 and ir/6 correspond to the three examples above. 
Let us mention that the class of isoradial graphs has been extensively studied in order to 
understand universality; see [21 [7] . 

Example 5. Let us conclude with one last example which does not belong to the class of 
isoradial weighted graphs. Let G denote the 2x1 square lattice with arbitrary coupling 
constants J\, J2, J3, J4 as illustrated in the right-hand side of Figure [TJ Then, the critical 
inverse temperature is given by the equation 

1 + £3X4 = X 3 + X4 + XlX 2 + XlX 2 X 3 + X2X 3 X4 + X1X2X3X4 

where Xi = tanh(/3Jj) for i = 1, 2, 3, 4. 

Organization of the article. The next section defines the Kac-Ward matrices and recalls 
several of their properties. Section 3 presents a connection between Kac-Ward and Kaste- 
leyn matrices, as well as consequences for the Ising model. Section 4 contains the proofs of 
Theorem 11.11 and Corollary 11.21 

Acknowledgments. The second author was supported by the ANR grant BLAN06-3-134462, 
the ERC AG CONFRA, as well as by the Swiss FNS. The authors would like to thank Yvan 
Velenik for very interesting comments and references. The authors also thank Zhongyang Li, 
Beatrice de Tilliere and Cedric Boutillier for helpful discussions. 
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Figure 3. The angle a(e,e') G (— 7r,7r). 



2. The Kac-Ward matrices 



The aim of this section is to review the definition and main properties of the Kac-Ward 
matrices associated with graphs embedded in surfaces. To simplify the present exposition, 
we shall only treat the case of toric graphs with straight edges, referring the reader to [9] for 
the general case and the proofs. 

Let us start with some general terminology and notation. Given a weighted graph (G, x), 
let E = E(G) be the set of oriented edges of G. Following [33], we shall denote by o(e) the 
origin of an oriented edge e G E, by t(e) its terminus, and by e the same edge with the 
opposite orientation. By abuse of notation, we shall write x e = x s for the weight associated 
to the unoriented edge corresponding to e and e. 

Now, assume that G is embedded in the oriented torus T 2 so that its edges are straight 
lines and its faces are topological discs. Fix a character ip of the fundamental group of T 2 , 
that is, an element of 



Note that the choice of two oriented simple closed curves jx,Jy representing a basis of 
-ffi(T 2 ;Z) determines an isomorphism ff 1 (T 2 ;C*) ~ (C*) 2 . Then, such a character simply 
corresponds to a pair of non-zero complex numbers (z, w). Let us now assume that the curves 
7x,1y avoid the vertices of G and intersect its edges transversally. Then, a natural 1-cocycle 
representing the class (z,w) is given by the map ip: E — > C* defined by </?(e) = z lx ' e w" /y ' e , 
where • denotes the intersection form. In words, 7 • e gathers a +1 (resp. —1) each time 7 
and e intersect in such a way that (7, e) define the positive (resp. negative) orientation on 



where a(e,e') G (— tt,tt) denotes the angle from e to e' , as illustrated in Figure O We shall 
call the matrix / — T v the Kac-Ward matrix associated to the weighted graph (G, x), and 
denote its determinant by P' p (G,x). When ip is identified with a pair of complex numbers 
(z,w), we shall simply denote the determinant by P z ' w (G,x). 



^Recall that in the present context, a 1-cocham is a map ip: E — !> C* such that <fi(e) = <^(e) _1 for all e G E. 
Such a map is called a 1-cocycle if for each face / of G C T 2 , <p(df) := Flees/ < < 5 ( e ) ~ Multiplying each 
93(e) such that o(e) = v by a fixed A„ G C* results in another 1-cocycle, which is said to be cohomologous to 
ip. Equivalence classes of 1-cocycles define the first cohomology group i/ 1 (T 2 ;C*). 



Hom(^!(T 2 ),C*) = H 1 ^ 2 -^) 
the first cohomology group of T 2 with coefficients in C* J3 



T 2 . 
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Expanding this determinant, see [9], leads to an expression of the form 

(3) F"(G,x) = Yl (-l) e(7) *(7M7), 

7er(G) 

where 6(7) € {0, 1} is some sign and elements of T(G) are unions of closed paths in G that 
never backtrack and pass through each edge of G at most twice, and if so, in opposite direc- 
tions. This shows that P lfi (G,x) does not depend on the choice of the 1-cocycle representing 
the cohomology class ip G H l (Y 2 ; £.*). Since 7 belongs to T(G) if and only if 7 does, we also 
immediately obtain that if ip belongs to H l (T 2 ] S l ) ~ S 1 x S 1 and x to WL E{ - G \ then P ip (G, x) 
is a real number. 

As detailed in [9], a closer analysis of the sign 6(7) leads to the following results. Fix two 
oriented simple closed curves j x , j y as above, thus identifying if 1 (T 2 ; C* ) with (C*) 2 and 
iii(T 2 ;Z2) with (Z2) 2 . For a S Hi(T 2 ; Z2) , let Z a denote the corresponding partial partition 
function, that is 

Z<x = 

where the sum ranges over all paths with homology class equal to a. 

Proposition 2.1 ([9]). We have the following equalities: 

P 1 ' l (G,x) = (Z 00 ~ Z w - Zoi - Z u ) 2 , P-^Cx) = (Zoo + Z 10 - Z 01 + Z n ) 2 , 

P^- x (G,x) = (Zoo ~ Z w + Zoi + Z u ) 2 , p-^iCx) = (Zoo + Z 10 + Z 01 - Z u ) 2 . 

In particular, this shows that for z, w = ±1, P Z ' W (G, x) is the square of a polynomial in the 
weight variables x. By construction, the constant coefficient (with respect to the variables x) 
of the polynomial P z ' w (G,x) is equal to 1. Denoting by P Z ' W (G, x) 1 / 2 the square root of the 
polynomial with constant coefficient +1, we get the following formula. 

Theorem 2.2 (|9j). The Ising partition function for a toric weighted graph (G,x) is given 
by 

Z (G, x) = 1 ( - P^iG, x) 1 ' 2 + P 1 '- 1 ^, xfl 2 + P-^(G, x) 1 ' 2 + P-^~\G, x) 1 / 2 ) . 

3. Relation to dimers and consequences 

The aim of this section is to show that the Kac-Ward determinant P Lp (G, x) is proportional 
to the Kasteleyn determinant of an associated bipartite weighted graph (Cc,y) (see Theo- 
rem [3J] below). By Kenyon and Okounkov's [214 Theorem 1], it follows that the curve defined 
by the zero-set of P' f (G,x) is a (simple) Harnack curve. It also implies some new avatar of 
Kramers- Wannier duality (Corollary 13. 3|) . 

Let (G,x) C T 2 be a weighted graph, and let us parametrize x e S (0, 1) by x e = tan(0 e /2) 
with 6 e E (0, 7r/2). Following Fan-Wu [2] and Dubedat [11], let us consider the associated 
weighted graph (Ca,y) C T 2 obtained from G as follows. Replace each edge e of G by a 
rectangle with the edges parallel to e having weight sin(# e ) while the other two edges have 
weight cos(# e ). In each corner of each face of G C T 2 , we now have two vertices; join them 
with an edge of weight 1. This is illustrated in FigureHl Note that since the torus is orientable, 
the graph Cq is bipartite: its vertices can be split into two sets B\JW (say, black and white 
vertices) such that no edge of Cq joins two vertices of the same group. 
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tan(0/2) 




C G 



sin(0) 




V 



sin(<9) 



' cos(0) cos(0) 
Figure 4. The weighted bipartite graph Cq associated to the weighted graph G. 



Recall that a Kasteleyn orientation \\.7\ [T8| [T9] on a bipartite graph C embedded in an 
orientable surface can be understood as a map oj: -E'(C) — >• {±1} such that for each face / of 
C, 



(&;)(/) := ]J w(e) = (-1) 



+ 1 



The associated Kasteleyn operator K(C, y) : C B — > C w can be defined by its matrix elements: 
for w G W and b € B, set 



= ^ w(e)y e , 



e:w— >b 

the sum being over all edges of C joining w to 6. If C C T 2 is endowed with a map 
p: E(C) —7- C* (for example, a 1-cocycle), then one can extend this definition to an op- 
erator K^iC^y): C B C w by multiplying the coefficient corresponding to the edge e by 
<p(e), where e is oriented from the white to the black vertex. 
The main result of this section is the following. 

Theorem 3.1. For any weighted graph (G,x) C T 2 , there is a Kasteleyn orientation on 
C G C T 2 such that 

P^(G,x) = 2-l y ( G )l \{ (1 + x 2 ) det(^(C G ,y)) 

eeE(G) 



for all ^ G iJ 1 (T 2 ; C) ~ (C*) 2 . 



Proof. The strategy of the proof is to gradually transform the Kac-Ward matrix for (G, x) 
into the Kasteleyn matrix K v (CG,y) while keeping track, at each step, of the effect of the 
transformation on the determinant. Let us begin with some notation. We shall write ££ (E) 
for the complex vector space spanned by the set E of oriented edges of G. Obviously, E 
can be partitioned into E = U«ev(G) wnere E v contains all oriented edges e with origin 
o(e) = v. Now, let us cyclically order the elements of E v by turning counterclockwise around 
v. (As T 2 is orientable, this can be done in a consistent way.) Given e G E v , let R(e) denote 
the next edge with respect to this cyclic order, as illustrated in Figure [5j This induces an 
automorphism R of «Sf(E). Also, let J denote the automorphism of Jz?(E) given by J(e) = e. 
Finally, we shall write x for the automorphism of Jz? (E) given by x(e) = x e e, and similarly 
for any weight system and for ip. 
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Figure 5. The endomorphism R of Jz?(E), and the angle j3 e . 

Let Succ G End(Jf(M.)) be defined as follows: if e is an oriented edge with terminus t(e) = v, 
then 

Succ(e) = (p(e)x e u(e,e')e, 

e'GE v 

where ui(e, e') = exp(|a(e, e')) for e' / e £ E„ with a(e, e') as in FigureEl and u)(e, e) = —i. 
Also, let T G End{££ (E)) be the endomorphism given by T = Succ + iJipx. By definition, 
the Kac-Ward matrix is equal to I — T. Now, consider the matrix 

A = (I — T)(J + iJtpx) = I — Succ + Com, 

where 

Com(e) = —icp(e)x e T(e) = —ix\ ui(e,e')e' 
if e has origin o(e) = u. Since 

det ( / + ^) = ndet ( i *® x < ) = n(i+*&, 

we get the equality 

(4) pv(G, x) = det (J - T) = JJ (1 + x^)" 1 det A. 

We now focus on the computation of det A. We shall transform the matrix A by multiplying 
it to the left with some well-chosen matrix N, whose determinant is easy to compute. Let 
Q G End(^(E)) be defined by 

Q(e) = exp (|/3 e ) e, 

where (3 e = tt — a(J(R(e)), e) G (0, 27r) denotes the angle between e and i?(e) (see Figure [5]). 
Obviously, the endomorphism N := I — i?C/ decomposes into iV = ©„ e y(G) iV w with iV^ G 
End{J£ \E V )) , and one easily computes 

detJVe = 1 - exp (§/3 e ) = 2. 
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Hence, the determinant of N is 2^ V ^ G ^. Now, let us compute the composition NA. If e has 
terminus t(e) = v, then 

NSucc(e) = (I — RQ)(p(e)x e o;(e, e')e' 

e'€E v 

= ip{e)x e ^2 M e > e ') - w(e,i?~ 1 (e'))exp (§/3.R-i( e '))) e' 

e'£E v 

= —2iip(e)x e e. 

Therefore, we have the equality NSucc = —2iJipx. Similarly, given e with origin o(e) = v, 
NCom(e) =(/ - RQ){-i)x 2 e ^ co(e,e')e' 

e'£E v \{e} 

= ~ ix l Yl ( w (e,e')-w(e, J R- 1 (e'))exp(§/3 iJ -i (e0 ))e' 

e>£E v \{e,R(e)} 

- ixl (^(e, R(ej) R(e) - u(e, # _1 (e)) exp (§/?fl-i( e )) e) 
= -xl(I + RQ)(e). 

These two equalities lead to 

NA = N(I - Succ + Com) 

= (J — RQ) + 2iJipx -(1 + RQ)x 2 
= (1 - x 2 ) + 2iJipx - RQ(1 + x 2 ). 

Since the determinant of N is 2^ V ^ G ^, the equality displayed above together with Equation Q 
give 

(5) P*{G,x) = 2~\ y ^\ [J (l + ^e)detM, 

e€E(G) 

where M is given by 

1 — x 2 2x 

M = = + iJip k - RQ = cos(0) + iJipsmiO) - RQ, 

1 + x z 1 + x z 

using the paramatrization x e = tan(0 e /2) of the weights. 

The final step is now to show that this matrix M is conjugate to the Kasteleyn matrix 
K^(Cg, y)- The theorem will then follow from Equation ([5]). Let ips '■ E — > B (resp. ipw : E — > 
W) denote the bijection mapping each oriented edge e of G to the unique black (resp. white) 
vertex of Cq immediately to the right (resp. left) of e, as illustrated below. Observe that the 
three maps i/jb ° > ° J ° and ips ° R ipyy associate to a fixed w € W the three 
black vertices of Cq adjacent to w. 



R(e) \i> B (R(e)) 



i> B {J{e)) 
• 



• ^s(e) e = W 1 H 
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Therefore, the operator 

= (ip B o M o t/v 1 )* : C B 
is given by the coefficients 

cos(0 e ) if the edge (w, b) is perpendicular to e 6 E; 

(p(e)isia(9 e ) if (iw, 6) is to the left of e £ E; 

— exp (|/3 e ) if (w, 6) is in the "corner" of e and i?(e); 

if w and 6 are not adjacent in Cq- 
This is illustrated below. 



R(e) y 




-exp(|/? e ) 
y(e)i sin(# e ) 



/* (p(e)- 1 zsin(0 e ) \ cos(6 e ) 



In other words, if^ is precisely the Kasteleyn operator of (Cg,u) associated to the 1-cocycle 
<p: E(C G ) -> C* given by 



<f(w,b) 



ip(e) if 6) runs parallel to e € E; 
1 else, 



and to the map Co: E(Cg) S 1 C C* given by 

{1 if (w, b) is perpendicular to e € E; 

z if (iw, 6) is to the left of e G E; 

— exp (§/3 e ) if (w, &) is in the corner of e and R(e). 

Since the 1-cocycles ip and induce the same class in ff 1 (T 2 ;C*), it only remains to handle 
the map Co. Extend it to a 1-cochain Co: E(C G ) — >■ S 1 by setting Co(b,w) = Co(w,b)~ l . Now, 
observe that for any face / of C G C T 2 , 

:= [J w(e) = 

This is obvious for the rectangular faces; for the faces corresponding to vertices of G, use the 
fact that the angles (3 e add up to 2ir around each vertex; for the faces corresponding to faces 
of G, use the fact that the angles a(e, e') add up to 2tt around each face. Furthermore, since 
we are working on the flat torus, one easily checks that Co(j) = ±1 if 7 denotes a 1-cycle in 
C G C T 2 . Therefore, Co is cohomologous to a Kasteleyn orientation 10. In other words, Co 
can be transformed into w by a sequence of the following transformation: multiply all the 
edges adjacent to a fixed vertex of C G by some complex number of modulus 1. Therefore, a 
Kasteleyn matrix i^(C G ,y) can be obtained from by multiplying lines and columns by 
complex numbers of modulus 1, and the equality 

-\V(G)\ 



P*{G,x) = T 



n <» 

e£E{G) 



det(K^(C G ,y)) 
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holds up to multiplication by a fixed complex number of modulus 1. For ip taking values in 
{±1}, both sides are real; therefore, the equality holds up to sign. One can then choose a 
Kasteleyn orientation such that the identity holds. □ 

Let us now fix a geometric basis of Hi(T 2 ;Z). As explained in Section [21 this allows to 
identify any element tp of i/ 1 (T 2 ;C*) with a pair of non-zero complex numbers (z,w) and 
to write P' p (G,x) = P Z,W (G, x). The Newton polygon of such a polynomial P z,w (G,x) = 
^2nmez a nmZ n w m is the convex hull of {(n,m) G Z 2 ; a nm ^ 0}. Recall that a weighted 
graph (G, x) C T 2 is non- degenerate if all the weights are positive and the complement of G 
in the torus consists in topological discs. Using (j3J), one easily checks that for such graphs, 
the Newton polygon of P z,w (G,x) has positive area. Theorem 13.11 and |21l Theorem 1] then 
immediately imply that (the real part of) the associated spectral curve 

A = {(z,w) G (C*) 2 ; P z ' w {G,x) = 0} C (C*) 2 

is a (possibly singular) Harnack curve. This means that the curve A C (C*) 2 intersects each 
torus T 2 (r, s) = {(z, w) G (C*) 2 ; \z\ = r, \w\ = s} in at most two points (see [30]). 

Recall that an element (zq,Wq) of a complex curve A = {(z,w) G (C*) 2 ; P(z,w) = 0} is 
called a singularity of A if J^P(z, w) = Jj^P(z, w) = 0. As shown in [3Q|, Lemma 6], Harnack 
curves only admit (a very specific type of) real singularities, meaning that both coordinates 
are real. In our case, this translates into the following corollary: 

Corollary 3.2. The spectral curve associated to a non- degenerate weighted graph embedded 
in the torus has no singularities other than real ones. 

We now turn to Kramers- Wannier type duality for the Kac-Ward determinants. Again, 
these results hold for the more general case of an arbitrary weighted graph embedded in a 
closed orientable surface. For the simplicity of this exposition, we shall only consider the 
special case of the torus. 

If G is embedded in the torus, its dual is the graph G* C T 2 obtained as follows: each face 
of G C T 2 defines a vertex of G* , and each edge of G bounding two faces of G C T 2 defines 
an edge between the two corresponding vertices of G* . Note that (G*)* = G. Finally, if G is 
endowed with weights x = (x e ) G (0, 1) E ( G \ define the dual weights x* = (x*) G (0, 1) E ( G ) via 
the condition x + x* + xx* = 1. If we use the parametrization x = tan(#/2), X* = tan(#*/2), 
then 6 and 9* are simply related by 9 + 8* = ir/2. Therefore, the weighted graph {Cq* , y(%*)) 
associated to (G*,x*) is equal to the weighted graph (Ca,y(x)) associated to (G,x). Hence, 

Theorem 13.11 together with the equality = immediately lead to the following. 

Corollary 3.3. For any toric weighted graph (G,x) and any tp G if 1 (T 2 ;C*) ; 

2 ' y(G) ' I I {t + XeT l P v {G,x) = 2l^ G *)l Y[ (l + ^) _1 -P v (G*,x*). □ 

eg£(G) eeE(G) 

As mentioned in Proposition 12.11 for (p = (z,w) G {±1} 2 , P lp {G,x) is the square of a 
polynomial in the weight variables x e . As the constant coefficient of P <p (G,x) is equal to 1, 
we can pick such a square root pv(G, x) 1 / 2 by requiring its constant coefficient to be +1. 
Taking a closer look at the sign leads to the following duality. 

Corollary 3.4. For any toric weighted graph (G, x) and any (p G -ff^CT 2 ; {±1}), 

2 \v(G)\/2 Yl (l + x e )" 1 / 2 P^(G,x) 1 / 2 = (-l) A ^2l v/ ( G *)l/ 2 J] (l + xtr^P^G^x*) 1 / 2 , 

e£E(G) e&E(G) 
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where A(<p) = 1 if tp = (1, 1), and A(cp) = else. 

Proof. By Corollary 13.31 we only need to determine the sign A{ip) in the equation above. 
Setting x = 1 (and therefore, x* = 0) leads to 

P^(G, 1) 1/2 = (-i)M<e) 2 (\ v ( G *)\+\E(G)\-\V(G)\)/2 = (_ 1 )A( V ) 2 |V(G*)| ) 

using the fact that |F(G)| — |i?(G)| + \ V {G*)\ is equal to the Euler characteristic of the torus, 
i.e. zero. Furthermore, the Ising partition function Z(G, x) with weights x = 1 is nothing but 
the cardinality of the Z2-vector space of 1-cycles modulo 2 in G. Since G is connected, the 
dimension of this space is classically equal to |-E(G)| — | V(G)| + 1 = | V(G*)| + 1. Theorem 12.21 
now reads 

2 |v(g-)|+i = i( - (-i)^ 1 ' 1 ) + (_i)^(L-i) + + (_i)^(-i -i))2l y ( G *)l. 

The term in parentheses is therefore equal to 4, a fact which determines the sign of the four 
terms. The corollary follows. □ 



4. The critical temperature 

Consider a planar non-degenerate locally-finite weighted graph (£f , J) invariant under a 
lattice A ~ Z Z. Alternatively, we will use (£f , x) when working directly with the high- 
temperature expansion. Recall that G = @/A. For integral positive n, m, let A nm ~ nL®mL 
and let G nm denote the toric weighted graph given by & /A nm . Note that G\\ = G. 

Let us introduce the free energy per fundamental domain as follows: 

(6) \ogZ x := lim log Z(G nn , x). 

n— yoo n 

This definition is justified by classical super-multiplicative properties of partition functions. 
The strategy of the proof of Theorem 11.11 is the following. We express the free energy in 
terms of the Kac-Ward determinant, and we show that for any x such that P z ' w (G, x) has 
a zero on T 2 := {(z,w) : \z\ = 1, \ w\ = 1}, the free energy \ogZ x is not twice differentiable 
at x. We then harness standard arguments on the Ising model to show that log Z x is twice 
differentiable except possibly at criticality. We begin by a classical lemma. 

Lemma 4.1. For any (x,w) € T 2 and x G R E ( G \ 

P z ' w (G nm ,x)= H H P u > v (G,x). 

u n =z v m =w 

Proof. The proof of [22, Theorem 3.3] applies almost verbatim. □ 

The next lemma shows that the only zeros of P z ' w (G,x) are localized at (1, 1). 
Lemma 4.2. For any (z, w) e T 2 \ {(1, 1)} and x G (0, l) E ^ G \ P Z ' W (G, x) > 0. 
Proof. The proof is divided into four steps. 

Step 1. P*> w (G,x) > for any (z,w) G {(-1, 1), (1, -1), (-1, -1)}. Proposition O shows 
that 

Z 00 (G, x) = i (P 1 - 1 ^, x) 1 ' 2 + PL-^G, x) 1 / 2 + p-^iG, x)V2 + p-i.-i(G, x) 1 ' 2 ) . 
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Corollary 13.41 and Theorem 12.21 then give the equality Zqq(G, x) = C ■ Z(G*,x*), where 

Q . = 2 |V(G*)|/2-|V(G)|/2-l J| fL 



1 + ^ 1/2 



+ X* 



In the same way, Z\q(G,x) can be expressed as a linear combination of P 1 ' 1 (G*, x*) 1 / 2 , 
P 1 '- 1 {G*,x*) 1 / 2 , P~ 1 ' 1 {G*,x*) 1 / 2 and P~ l ~ l (G*,x*) 1 / 2 . Using Proposition EH again, we 
obtain the equality 

Z 10 (G,x) = C-(Z 00 (G*,x*) + Z 10 (G*,x*)-Z 01 (G*,x*)-Z n {G*,x*)), 

which leads to Z\q(G,x) < C ■ Z(G*,x*) = Zqq(G,x). This argument can be carried out for 
any homology class a, so 

Z a (G,x) < Z 00 {G,x) for any a G {00,01,10,11}. 

The assumption that G is non-degenerate implies that all Z a (G, x)'s are strictly positive. The 
statement now follows from the inequality displayed above, and Proposition 12. H 

Step 2. P Z ' W (G, x) / for any (z, w) such that z m = — 1 and w n = — 1 for some (m, n) G N 2 . 
This follows immediately from the first step applied to G mn and Lemma 14.11 

Step 3. P Z ' W (G, x) > for any (z, w) G T 2 . Recall that P Z ' W (G, x) is real for any (z, w) G T 2 . 
Let us fix (z,w) with z n = —1 and w m = —1. By the second step, 

X = {x G (0,l) E(G) I P z ' w (G,x) > 0} = {x G (0,1) B ( G ) | P z ' w (G,x) > 0}. 

By continuity of x i— > P z,w (G,x), X is therefore both closed and open. It is also non-empty 
since P Z ' W (G, x) tends to 1 as x tends to 0. By connexity of (0, 1) E ( G \ X is this whole space. 
By continuity of (z,w) \— > P z ' w (G,x), the statement follows. 

Step 4. P z > w (G,x) = implies (z,w) = (1,1). Assume that P z > w {G,x) = 0. Then, (z,w) 
must be a singularity (i.e. satisfy J^P z ' w (G,x) = J^P z ' w (G,x) = 0); otherwise, P z ' w (G,x) 
would take negative values near (z, w), contradicting the third step. By Corollarv l3.2l (z, w) G 
T 2 is real, i.e. (z,w) G {(— 1, — 1), (— 1, 1), (1, — 1), (1, 1)}. By the first step, (z,w) must be 
equal to (1,1), and the lemma is proved. □ 

The free energy can be expressed in terms of P as follows: 

Lemma 4.3. For any x G (0, l) E<yG \ 

l ogZx = — 1— [ logP z > w (G,x) — — . 
2(2iri) z Jj2 z w 

Proof. The proof is inspired by the proof of |22i Theorem 3.5]. First note that P Z ' W (G, x) > 
for any (z,w) ^ (1, 1) by Lemma 14.21 an d that 

logP z ' w (G,x) = 0[log(|* - 1| + \w- 1|)], 

which legitimates the integral on the right-hand side. Lemma 14.11 and the bounded conver- 
gence theorem imply that for (e, 77) G {(—1, —1), (—1, 1), (1, —1)}, 

±logP^(G nn ,x) = A E E logJ^' w (G,x) — > -^-2 / \ogP^(G,x)- — . 

zn=EW n =v v 1 J 1 

Now, Proposition 12.11 and Theorem 12.21 imply the inequalities 

p- x > x (G nn ,x) < Z(G rm ,x) 2 < Q/4m a x{P- 1 '- 1 (G nn ,x),p- 1 ' 1 {G nn ,x),P 1 '- 1 {G nn ,x)}, 
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which lead to the claim. □ 

For the next theorem, let us adopt the terminology of the high-temperature expansion. For 
(&,x) biperiodic, let fj,^ be the Ising measure on & with edge-weights x and + boundary 
conditions. Let x* such that x + x* + xx* = 1 be the dual weights obtained by Kramers- 
Wannier duality. 

Theorem 4.4. Let ^ be a non- degenerate locally-finite doubly periodic graph and r G V(SP). 
Then, 

(i) If [i^ y (a r ) = for any weights y in a neighborhood of x, then there exists c = c(x) > 
such that \im x {p a Gb) < exp(— c\a — b\), for any a, b G V(SP); 

(ii) If fJ>@ y (o- r ) > for any weights y in a neighborhood of x, there exists c' = c'{x) > 
such that fig* x *(o~ u a v ) < exp(— c'\u — v\), for any u,v G V(@*). 

While this theorem is not surprising and follows from very classical ingredients, the proof 
does not appear in the literature. We therefore recall it here. 



Proof of Theorem \4-4\ Let us prove (i). Choose (3 and J in such a way that x e = tanh(/3J e ). 
The condition implies that j3 < (3 C for (3f,J). We conclude by harnessing pQ Theorem 1]. 
Note that this theorem applies in the very general context of finite range Ising models on "L d 
with periodic coupling constants. In our case, the model is only biperiodic but as discussed 
by the authors, the proof extends very easily to this framework. 

Let us now deal with (ii). We aim to apply (i) to the dual measures. For this rea- 
son, it is sufficient to prove that /xj y { a r) > implies y *{o- u ) = or equivalently that 
Liyt e y *(o~ u cr v ) — > as \u — v \ — > oo (u, v G Sf*), where "free" refers to free boundary conditions. 
Intuitively, this claim is valid since there cannot be a positive spontaneous magnetization for 
both the primal and dual Ising models. This is best seen in the context of random-cluster 
models. We thus harness the Edwards-Sokal coupling |13j . 

Let <py 2 De the random-cluster measure on <$ with cluster-weight 2, edge-weights given 

by 

1 + Ve 

and wired boundary conditions; see [HI Section 4.2]. The Edwards-Sokal coupling [161 Section 
1.4] shows that 

4>h,p,2( r i s connected to infinity) = pi^ y {o~ r ) > 
which implies the existence of an infinite cluster (fy 2 -almost surely. Let (<f>^ 2 )* be the dual 
measure, see |16|. Section 6.1]. Since the random-cluster model satisfies the FKG inequal- 
ity [161 Theorem 3.8], Corollary 9.4.6 of [34] implies that there cannot be coexistence of an 
infinite cluster and an infinite dual cluster 0^ 2 -almost surely. Thus, there is no infinite clus- 
ter on <$* for the dual random-cluster model (4>^ 2 )*-almost surely. Now, (cj>y 2 )* = ^» e p * 2 
and Hcffy* are also coupled via the Edwards-Sokal coupling. We obtain 

fj,^T y *(a u a v ) = (0^ iP)2 )*(^ connected to v) — > 
as | it — v\ — > oo (n, v G □ 

Remark 4.5. For (ii), one can also invoke (with some modifications) the result of Lebowitz 
and Pfister [23] together with duality. 
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Proof of Theorem \l.l[ Fix J). Define x$ = (tanh(/3J e )) e where j3 > 0. First, Proposi- 
tion O implies that Equation © is equivalent to P 1 ' 1 (G, xp) 1 / 2 = 0. Note that P 1 ' l (G, xp) 1 / 2 
tends to 1 as f3 tends to (by definition), and to —2\ V ( G *'\ as j3 tends to oo (by Corollary 13. 4|) . 
so that there exists at least one solution (in j3) to this equation. We now show that there 
exists a unique such solution by proving that P 1 ' 1 (G, x@) = implies f3 = j3 c . 

We first show that (5 C < (5 by assuming that /3 < /3 C , or equivalently x^ < x^, and by seeking 
for a contradiction. Since P 1 ' 1 (G, x) is not constant on a neighborhood of x# (for instance 
P 1 ' 1 (G, x^/) > for p' close enough to /3), there exist e G -E(G) and x^ < x < x^ c such 
that P 1,1 (G, x) = and x e i— )■ P 1 ' 1 ^, x) is non constant. Fix such an edge e and weights 
x, and let x(t) be defined by x(f) e / = x e / if e' 7^ e and x(i) e = x e + t. By Equation ([3]), 
t I—)- P 1,:L (G, x(t)) is a polynomial of degree exactly 2. Furthermore, since P 1,1 (G, x(0)) = 0, 
pe 1 ,e (G^x(t)) = P e ' e 1V (G, x(t)) and P 6 ' ,e (G,x) > for any (9,rj,t) in a neighborhood 
of the origin (Lemma 14. 2p . we obtain the following development near (0,0,0): 

P ei6 ' eiv (G,x(t)) = (a u t 2 + a 22 6 2 + a 33V 2 + a 23 6 V ) f (t , 6 , 77) , 

where f{t,6,rj) = 1 + o(|i, #,fy| 2 ) is a non-vanishing analytic function, and the coefficients 
satisfy 022033 — 3023 — and an > 0. Now, 



log/(t,M)d0d77 

is twice differentiable in £, so that \ogZ x is twice-differentiable in x e at x if and only if 

/7T /*7T 
/ log(ant 2 + a22^ 2 + a33^? 2 + a 23 8rj)d9dr) 
-TV J —IX 

is twice differentiable at 0. For t 7^ 0, the second derivative in t of this function equals 

2a 11 (a 11 t 2 + a 22 6 2 + a 33 r] 2 + a 23 6r]) - \a\ x t 2 ^ 



(ant 2 + a 22 9 2 + a 33 r] 2 + a 23 9ij) 
As t tends to 0, this integral tends to 

2an 



-7T /-tt (a22^ 2 + 033?? 2 + 023^??) 



dtfdr/, 



that is, to 00, since an > and 022^33 — 3 a 23 — 0- But this is in contradiction with the 
assumption that x < x^ c since in this case, exponential decay implies that log Z x is twice 
differentiable. We now justify this last statement. Fix a representative {a, b} £ E(@) of the 
edge e. Let & nn be a fundamental domain of the action of A nn on <S . Further assume that ^ nn 
contains the edge e. Let E n (resp. E) be the set of translates of e in E(& nn ) (resp. E(&)). 
Since the definition of the free energy does not depend on the boundary condition, one has 

logZ x = lim —^\ogZ(% n ,x), 

n— >oo n 

where Z(& nn ,x) is the partition function on ^Snn with, free boundary conditions. Set x e — 
tanh(/3Jg). The high temperature expansion ([T]) shows that 

logZ, = lim -llogZ/'(Sf nn ) - V log(cosh(/3^)) - |V(G)| • log(2), 

eg£(G) 
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where Zi ' (@ nn ) is defined in the introduction. Since J' e depends smoothly on x e , it is sufficient 
to show that liuin^oo -^logZl (& nn ) is twice differentiable with respect to J' e . We obtain 



J /3 



2 dJ! 2 



71- 



{u,v}eE n 



where ^ nn , x is the measure on with free boundary conditions. Lebowitz's inequality 
Remark (i), p. 91] then yields 

The comparison between boundary conditions implies lA<g nn A a aO~u) < x { a a^u), which, 
together with Lebowitz's inequality and Property (i) of Theorem 14.41 leads to 

It^nnA^bVuPv) ~ (J# nn A <7 a (T b)lMfn n A (T u<Tv) <2 ^ exp(-c(a;) \u - a\ ) . 

{u,v}eE n {u,v}GE 

The term on the right-hand side is therefore bounded uniformly in n and in x, provided that 
x takes value in a compact subset of {x < X/? c }. The bounded convergence theorem then 
implies that \ogZ x is twice differentiable in J' e , and thus in x e . In conclusion, x cannot be 
smaller than xa, and we obtain that (3 > /3 C . 



Let us conclude the proof by showing that /3 < /3 C . Corollary [373] shows that P 1 ' 1 (G, x^) = 
implies = 0. Now if j3 > (3 C , or equivalently x*^ < x*^ c , one can run the previous 

argument for the dual model, using Property (ii) of Theorem 14.41 in place of Property (i). □ 

Proof of Corollary [TE When f3 / /3 C , P z ^ w {G,xp) > for any (z,w) G T 2 and /3 ^ 
log P Z,W (G, xr) is analytic. By Lemma [4.31 the free energy is a parameter-dependent integral 
which is analytic at ft ^ f3 c . □ 
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